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Relativistic and non-relativistic fluid equations can exhibit finite time singular solutions including 
density singularities appearing in collapse or compression systems and gradient singularities in shock 
waves. However, only the non-relativistic fluid equations have been shown to exhibit finite-time 
velocity singularities, u — > oo in this regime. As this limit violates the cosmic speed limit v < c, 
where c is the speed of light, it is unclear how velocity singularity source terms in the non-relativistic 
equations are modified in the relativistic equations. This paper uses a driven system to demonstrate 
that singular points within the relativistic fluid equations can generate finite time relativistic velocity 
singularities i; — > c. 



I. INTRODUCTION 

The relativistic and non-relativistic fluid equations 
must necessarily provide identical descriptions of fluid 
motion in the non-relativistic regime suggesting that 
density and gradient singularities present in the non- 
relativistic regime must also be appear in the relativistic 
equations. In particular, density singularities p — > oo fea- 
ture in the relativistic gravitationally collapse of a non- 
rotating zero pressure dust cloud [Q, as well as in non- 
relativistic compressive systems such as shock tunnels 
where growing pressures burst thick metal diaphragms 
1^, in the hand-held Diamond- Anvil generating 2 mil- 
lion atmospheres of pressure , in cavitation generating 
energy densities sufficient to pit steel [|| , in sonolumines- 
cence generating temperatures in excess of 100,000 de- 
grees 1^, ^ 01 , in explosive compression of magnetic fields 
generating fields up to 1,000T || and in supernovae 
collapse experiments where inertial confinement gener- 
ates temperatures in excess of 100 million degrees 
Of course, in the non-relativistic regime fluid parameters 
do not literally become singular — the approximations de- 
riving the non-relativistic equations break down long be- 
fore that point. 

Similarly, gradient singularities appear within both the 
relativistic and non-relativistic fluid equations to gener- 
ate shock waves within relativistic jets |ll|, p^ , in non- 
relativistic supersonic flows I, m |l|, [l§7 and interest- 
ingly, within trombones p6[ |. 

Singular points present within the fluid equations can 
potentially drive not only density and gradient singu- 
lar solutions, but also flnite time velocity singulari- 
ties. In the non-relativistic regime, velocity singulari- 
ties t; — > oo have been exhibited in fissioning fluid drops 
| p7| , [ig|j |2C |_^hin-films evaporating or pinched to zero 
thickness [pi[, 22|, and curvature collapse and jet eruption 
from fluid surfaces These systems typically feature 
a free fluid surface of microscopic scale in one dimen- 
sion and macroscopic scale in other dimensions. This 
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scale disparity allows surface tension to dominate other 
forces in the microscopic dimensions allowing self-similar 
solutions where essentially, fluid is evacuated across some 
macroscopic distance on infinitesimal timescales generat- 
ing a singular velocity solution. These systems can be 
considerably simplified and scaled up in size by replacing 
the free fluid surface by a flat moving boundary which 
forces fluid evacuation over macroscopic distances on in- 
finitesimal timescales to again generate singular veloci- 
ties. The validity of this approach was experimentally 
confirmed in Ref. p^ . 

More generally, it is possible that similar singular 
source terms might support finite time current singu- 
larities in magnetohydrodynamic models of magnetic 
reconnection [ p5[ , and that unexplained car tire noise 
might result from nonlinearities present in such driven 
collapse systems Additionally, velocity singular- 

ities appear at boundary layer separation points [ p^ , 
while examinations of non-relativistic fluid singulari- 
ties and the hydrodynamic blow-up problem appear in 

13, H, 1^, 1^, |i|, g gg. 

To date, there has been no examination of the poten- 
tial for singular points within the relativistic equations 
to drive "relativistic velocity singularities" v —^ c, where 
here, c is the speed of light recently termed "Einstein's 
constant" Existing schemes to accelerate particles 

to near light speed do not exploit velocity singularities 
and typically rely on density singularities where essen- 
tially, "flreballs" accelerate particles as when coalescing 
neutron stars source gamma ray bursts ||36|| , asymmetric 
(or turbulent) supernovae collapse driving bullets of core 
material through overlying stellar layers p7| , |3^ , and high 
intensity table-top lasers drive particle acceleration [ p9[ . 
Other proposed mechanisms employ gradual particle ac- 
celeration over macroscopic distances and include shock 
wave surfing mechanisms generating cosmic ray particle 
energies up to lO^^eV and magnetohydrodynamic 
acceleration of polar jets from accretion disks 

In contrast to these approaches, this paper examines 
whether singular points within the relativistic fluid equa- 
tions can be exploited to generate near relativistic flows. 
The intractable nature of the problem requires excessive 
simpliflcation and, as in the non-relativistic case p4[|, we 
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exploit driven boundaries to force the relativistic fluid 
evacuation over macroscopic distances on infinitesimal 
timescales, and demonstrate in the appropriate limit that 
singular points within the relativistic fluid equations can 
drive fluid velocities to the limit v ^ c. 
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FIG. 1: (a) A schematic of a near relativistic column of per- 
fect fluid impacting an incompressible surface, and (b) mod- 
elling the column as a piston compressing underlying layers 
of gas. 



II. RELATIVISTIC VELOCITY 
SINGULARITIES 

Our initial examination of relativistic velocity singular- 
ities requires scaling up non-relativistic fluid evacuation 
systems to relativistic speeds. To this end, consider a 
perfect fluid column travelling at near relativistic speeds 
hitting an incompressible surface and generating an out- 
flow as shown in Fig. |l](a). The impacting fluid column 
acts as a compressive piston as shown in Fig. ^(b) so that 
as the piston closes, entrapped fluid must evacuate over 
macroscopic distances on ever shrinking timescales (as- 
suming the fluid remains relatively uncompressed.) This 
crude model suffices to allow the approximate solution of 
the fluid expulsion speed in horizontal directions. 

For simplicity we ignore gravity so fluid motion is 
entirely governed by the relevant conservation equa- 
tions, and consider a flat space-time with metric rjap = 



a/3 



diag (—1, 1, 1, 1) where primed coordinates x'' 



n 

(t' , x' ,y' , z') give the contravarient fluid velocity vector 
as U'" — (7,7v') with v' — {vx' ,Vy' ,Vz') and 7 — 
l/Vl — w'^ in geometrized units setting the speed of light 
to unity. (Greek indices run over values — 4 while Latin 
indices vary from 1 — 3.) 

A perfect gas has particle flux and energy-momentum 



N"^ = nU'" 



(1) 



where n, p and p are the momentarily comoving frame 
particle number density, pressure and energy density re- 
spectively Q. Conservation equations are N'^./^ = 
and T'f"".^ = giving 



D' (771) 

D'[j^{p + p)V] +V'p 



(2) 



with D' = df + V'.v' and V = (d^' , dy, , d,,). These 
equations give respectively, particle number conserva- 
tion, mass continuity and momentum conservation. 

The presence of a moving boundary (the descending 
piston) complicates analysis and it is useful to use time 
dependent computational space coordinates in which the 
descending piston is rendered stationary. If the piston is 
at height do at time < = 0, then its subsequent height 
can be written y'{t') = f{t')do in terms of a monotoni- 
cally decreasing function f{t') with /(O) = 1. Full closure 
occurs in the limit / — > 0, and the piston's vertical veloc- 
ity is Vy(t') = f{t')do. The transformation to unprimed 
computational space coordinates is then 



(3) 



so effectively, vertical compression in physical space 
Ay' corresponds to a constant separation in com- 
putational coordinates, Ax constant. The vertical com- 
pression is fully captured in the function f{t) = f{t'), 
and its effects can be assessed by tracking this function 
in governing equations. 

The new coordinates determine a new metric g^^ and 
proper time dr^ = g^^dx^dx'^ 




dr' = 



dt" + 2ffxdt dx 



Jx) 1 

+dx^ + fdx^ + dz^, 
while the non-zero affine connections are 



(4) 



(5) 



As required for a flat spacetime, the metric tensor g^i, has 
three positive and one negative eigenvalue and the Ricci 
tensor is zero everywhere. The transformed contravariant 
velocity vector is J7" = [7, 7 (u^;, fp^,, Wz)] with 



- fx 



(6) 
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Then, as required, the piston velocity is = as Vyi = 
fxi^ fdo). 

In the new coordinates, the particle flux and energy- 
momentum tensors become 



(7) 



SO the conservation equations N^^-^ = and T^";^ = 
give 

D{fin) = 
D{fj\p + p)v)+\/{fp) = { A \ , (8) 




with 



A 



fP 



P 



(9) 



and D ^ dt + V.v, V = {dx,d^,dz 



and 



l^iP + P) -P 



T*^ = 7^(P + + JXP- 



(10) 



These equations give respectively, particle number con- 
servation, mass continuity and momentum conservation 
as can be seen by taking the limit / = 1 and / = 0. 

An understanding of the physics embodied in Eqs. (|^) 
and (|^) can be obtained by considering solutions in two 
exactly solvable physical limits. The first considers a sta- 
tionary piston with the fluid in equilibrium. In this case, 
Eq. (g) is solved by setting fluid velocity v' = 0, and 
number density, pressure and density everywhere con- 
stant. Eq. reduces to Eq. (||) in the case f ~ 1, while 
for / < 1, Eq. (|) is solved by setting v = (0, -fx/f, 0), 
and again, uniform number density, pressure and den- 
sity. The second case considers a collisionless relativis- 
tic dust with p = (as some equation of state must 
be specified), and treats a closed piston lacking horizon- 
tal flows Vx' = Vz' ~ with the piston descending at 
a constant velocity (/ — 0). The dust at height y has 
initial vertical velocity at time i = of Vy' = fy and 
these dust elements preserve this velocity for all times as 
pressure is zero. For later times t at fixed height y, the 
fluid velocity increases as Vy' — fy/f- Eq. (^ is then 
solved by v' = (0, fy/f, 0) and everywhere uniformly in- 
creasing number density jnit) — ci/ f and fluid density 
j^p{t) = C2// for constants ci and C2. Conversely, the 
moving coordinate system tracks the initial dust elements 
leading to a constant velocity v = (0, 0, 0) which together 
with uniformly increasing number density "fn{t) — ci/ f 
and fluid density ^"^pit) = C2// solves Eq. (||). 



We now turn to consider the role of the squeezing pa- 
rameter f{t) in Eq. (^) which is linked everywhere on the 
left-hand-side (LHS) with 7. As is well known, the singu- 
larity 7 — > c» ensures that all fluid velocities remain less 
than the speed of light, v < c. Then, the countervailing 
linkage between this singularity and the zero asymptote 
of the squeezing parameter / — s- suggests a minimiza- 
tion of the effects of the combined terms /7 and /7^. 

The physical meaning of the squeezing parameter / 
can be seen in more detail using the injection form of the 
relativistic fluid equations 



D[^^{p + p)\-dtp 
D [7^(p + p)v] + Vp 



dxixp) 




{l^iP + P) -p) 



h\p + p) 

f 



(11) 



The left hand sides (LHS) here are identical to those 
of Eq. (||) while the right hand sides (RHS) represent 
sources injecting matter and momentum, including into 
the horizontal directions x and z. (Similar equations are 
used to describe rockets employing real injection systems 
to generate high speed expelled flows Q.) The injection 
sources are dominated in the relativistic large 7 regime 
by terms proportional to the ratio /// and either 7 or 7^. 
Then, as long as / ^ 0, the respective limits 7 — > 00 and 
/ ^ reinforce each other to inject (potentially) infinite 
mass and momentum in finite time to drive a horizontal 
relativistic velocity singularity. 

A very approximate analysis of the velocity and energy 
of the horizontal expelled jet can be performed but solu- 
tions are suggestive only. Consider the case where again 
the moving coordinate system approximately tracks the 
descending fluid so ~ 0, though now, the increasing 
pressure drives a near relativistic horizontal velocity Vx 
constituting an expelled jet. (Symmetry allows ignoring 
the z direction.) Fluid velocities are then (ua:,0,0). In 
the large 7 limit each of the above conservation equations 
gives the time rate of change of the expulsion velocity as 



dtVx 
dtVx 
dtVx 



V.V771 -t- "fh 

njHx 



f 



fl'^Vx 



v.WR + i'^R f{R~p) 



2-f^Rvx 
V.w-f'^Rvx 



2fRi^Vx 

dxP fVx 



72i?(l + 2^^vl) /(I + 272^2) 



(12) 



where R — ^'^{p + p) and the top line results from the 
conservation of particle number, the second line from 
the conservation of mass and the third line stems from 
the conservation of momentum in the x direction. The 
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asymptotic form of each of these equations in the large 7 
limit is approximately 



dtVx 



f 



(13) 



and this form is used for further discussion. The first 
term on the RHS is independent of / and, in the absence 
of squeezing (/ = 0), necessarily goes to zero so that 
dtVx ^ as 7 — > 00 ensuring v < c. The second term 
reflects squeezing effects and is positive (/ < 0) and the 
linkage between the limits / — > and 7 ^ 00 offsets the 
effect of 7. This equation can be integrated by ignoring 
the small first term with initial conditions /(to) = 1 a-nd 
Vxih) = Wo giving 



vl{t) = l-f{l-vl) 



or equivalently 



(14) 



(15) 



It is of interest to consider the non-physical mathematical 
limit / — ^ for some finite time T > to representing full 
closure so /(T) = 0. This hmit gives 



Vx{T) = 1, 



(16) 



a light speed jet. It is of course impossible for any ma- 
terial fluid to reach light speed though squeezed system 
might be able to exploit relativistic velocity singularity 
source terms to approach this ideal using singular mass 
and momentum injection source terms proportional to 

-hVf- 

Particle energy within the expelled jet is then 



E cx '-fit) (X y ^ 00, 
though this result is suggestive only. 



(17) 



Here, the diagonal pressure terms in the vertical direc- 
tion are scaled by a factor 1//^ over those in horizontal 
directions so that the limit / ^ creates the large pres- 
sure gradients driving the horizontal velocity singularity. 
Conservation of the energy- momentum tensor T^'^ = 
then gives the modified Euler equations 



dtU + dxF 



(19) 



fpVx 

f {pvl + p^h) 

fpVxVy 



fP 

U = I fpvx I F 

fpVy 
fP^X 

G = I fpVxVx 

fpVyV^+p'</-i 



Here = {vy — fx)l f a-^d we show mass continuity 
(top line) and momentum conservation in the x and y 
directions. (Symmetry constrains fluid motions to this 
plane.) Adiabatic perfect gases with p — p"' {'J = 1.4) 
are considered and all variables are dimensionless with 
dimensioned (primed) variables being given by x' = xL, 
v' = vuq, p' — ppq, p' — ppq and t' = tL/ao with L 
being some convenient length parameter and = 7Po/po 
giving the local speed of sound. 

Analytic solutions can be obtained in the low velocity 
limit Vx 1 (allowing v'^ « 0), and by assuming linear 
squeezing f{t) = 1 — t/T for some closure time T, with 
fluid elements possessing velocity Vy — fx tracked by the 
time dependent coordinates giving — 0. The analytic 
solution for an incompressible fluid with p — 1 (to give 
the maximum expulsion velocity) is then 



(x,t) 



fx 
/ 



(20) 



III. NON-RELATIVISTIC VELOCITY 
SINGULARITIES 

For completeness, details of the non-relativistic veloc- 
ity singularity created by asymmetrically compressing a 
fluid are given. (See Ref. for more details.) The 

conservative form of the inviscid and dimensionless Eu- 
ler equations with zero conductivity are derived from the 
the relativistic equations (0) after taking the small veloc- 
ity limit 7 — !■ 1, p ~ 0{pv'^, and p + p ^ p, so now p and 
p become the usual fluid density and pressure. The use 
of time-dependent coordinates [Eq. (||)] gives the Energy- 
Momentum tensor for a perfect gas as 

- p 

= pv^ 

= p«V-Hdiag(p,p//2,p). (18) 



displaying a velocity singularity as / — *■ 0. This solution 
is strictly valid only while Vx < 0.3 where an unconfined 
compressible fluid remains approximately uncompressed. 
Experiments demonstrate that this solution provides a 
reasonable approximation to the expulsion velocity into 
the supersonic regime [ p^ . 



IV. CONCLUSION 

This paper uses a preliminary and necessarily crude 
analysis to establish that relativistic velocity singulari- 
ties can appear in relativistic hydrodynamic equations. 
Velocity singularities provide such an elegant and simple 
means of sourcing near light speed relativistic jets that it 
would be surprising if this mechanism was not accessed 
in some astrophysical systems. 
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